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ABSTRACT 

The Palomar Transient Factory (PTF) is a fully-automated, wide-field survey aimed at a systematic explo- 
ration of the optical transient sky. The transient survey is performed using a new 8.1 square degree camera 
installed on the 48-inch Samuel Oschin telescope at Palomar Observatory; colors and light curves for detected 
transients are obtained with the automated Palomar 60-inch telescope. PTF uses eighty percent of the 1.2-m 
and fifty percent of the 1.5-m telescope time. With an exposure of 60-s the survey reaches a depth of m g /^21.3 
and mR^20.6 (5<x, median seeing). Four major experiments are planned for the five-year project: 1) a 5-day 
cadence supernova search; 2) a rapid transient search with cadences between 90 seconds and 1 day; 3) a search 
for eclipsing binaries and transiting planets in Orion; and 4) a 3n sr deep H-alpha survey. PTF provides auto- 
matic, realtime transient classification and followup, as well as a database including every source detected in 
each frame. This paper summarizes the PTF project, including several months of on- sky performance tests of 
the new survey camera, the observing plans and the data reduction strategy. We conclude by detailing the first 
51 PTF optical transient detections, found in commissioning data. 

Subject headings: Astronomical Instrumentation, Data Analysis and Techniques, Supernovae, Quasars and 
Active Galactic Nuclei, Stars 



1. INTRODUCTION 

The Palomar Transient Factory (PTF) is a comprehensive 
transient detection system including a wide-field survey cam- 
era, an automated realtime data reduction pipeline, a dedi- 
cated photometric follow up telescope, and a full archive of 
all detected sources. The survey camera achieved first light 
on 13 Dec 2008; the project is planned to complete commis- 
sioning in June 2009. 

This paper describes the technical aspects of the PTF 
project. An accompanying paper (Rau et al. 2009) describes 
the science planned for PTF in detail; subsequent papers will 
discuss the various PTF pipelines in more detail as well as the 
results from the PTF surveys. 

The transient detection survey component of PTF is per- 
formed at the automated Palomar Samuel Oschin 48-inch 
telescope (P48); candidate transients are photometrically 
followed up at the automated Palomar 60-inch telescope 
(P60). This dual-telescope approach allows both high survey 
throughput and a very flexible followup program. 

PTF's survey camera is based on the CFH12K mosaic cam- 
era formerly at the Canada-France-Hawaii Telescope, newly 
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mounted on the P48 telescope at Palomar Observatory, Cali- 
fornia (and called "the PTF Survey Camera"). The camera has 
101 megapixels, 1-arcsec sampling and a 8.1 square-degree 
field of view. Observations are mainly performed in one of 
two broad-band filters (Mould-R, SDSS-g'). Under median 
seeing conditions (1.1 arcsec at Palomar) the camera achieves 
2.0 arcsec FWHM images, and reaches 5<x limiting AB mag- 
nitudes of nig/^21.3 and niR^20.6 in 60-second exposures. 

Data taken with the camera are transferred to two auto- 
mated reduction pipelines (Figure Q]). A near-realtime im- 
age subtraction pipeline is run at Lawrence Berkeley National 
Laboratory (LBNL) and has the goal of identifying optical 
transients within minutes of images being taken. The output 
of this pipeline is sent to UC Berkeley where a source clas- 
sifier determines a set of probabilistic statements about the 
scientific classification of the transients based on all available 
time- series and context data. 

On few-day timescales the images are also ingested into 
a database at the Infrared Processing and Analysis Center 
(IPAC). Each incoming frame is calibrated and searched for 
objects, before the detections are merged into a database 
which is comprehensive, made public after an 1 8 month pro- 
prietary period and can be queried for all detections in an area 
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Fig. 1. — An overview of the PTF project data flow. 



of the sky. 

Followup of detected transients is a vital component of suc- 
cessful transient surveys. The P60 photometric followup tele- 
scope automatically generates colors and light curves for in- 
teresting transients detected using P48. The PTF collabora- 
tion also leverages a further 15 telescopes for photometric 
and spectroscopic followup. An automated system will col- 
late detections from the Berkeley classification engine, make 
them available to the various follow up facilities, coordinate 
the observations, and report on the results. 

This paper is organized as follows: in Section[2]we describe 
the new PTF survey camera and automated observing system, 
and detail the on-sky performance of the system. Section [3] 
discusses the initial PTF observing strategy. Section [4] de- 
scribes the automated data reduction pipeline and transient 
classification system, and Section \5\ details the PTF followup 
systems. In Section we conclude by describing PTF's first 
confirmed optical transient detections. 

2. THE PTF SURVEY SYSTEM: P48 + THE PTF SURVEY CAMERA 

The P48 telescope is a wide-field Schmidt with a 48-inch 
aperture, a glass c orrector plate and a 72-inch (f/2.5) mirror 
( Harrington! 1 1 952h . The telescop e performed both Palomar 
Sky Surveys (POSS I and POSS IIi lReid et al.l (ll99hY ). and re- 
centl y completed the Palom ar-Quest digital synoptic sky sur- 
vey dDjorgovski et al.l l2008). before the start of PTF modifi- 
cations in autumn 2008. 

The P TF survey camera (Fig ure|2]) is based on the CFH12K 
camera dCuillandre et al.ll2000h . The camera was extensively 
re-engineered by Caltech Optical Observatories for faster 
readout, closed-cycle thermal control, and robust survey op- 
eration. The CCD focal plane was untouched but the rest of 
the camera was modified for reliable and low-cost operation 
on the P48. Primary requirements for the design, in addition 
to the mechanical and optical modifications required for oper- 
ation on P48, were to reduce the operation cost of the camera, 
to minimize the beam obstruction, and to increase the readout 
speed for PTF operations. 

The P48 telescope optics are significantly faster than the 
CFHT, leading to stringent requirements on the camera's 
CCD array flatness and optical quality. The low-operation- 
cost requirement led to swapping the camera cooling system 
from a LN2 system to a CryoTiger closed-cycle cooler. The 
camera was also upgraded with a new precision shutter and 
filter-changer assembly. The telescope was refurbished for 
operation with PTF and a new queue- scheduling automated 
observatory control system was implemented. 



TABLE 1 

The specifications of the PTF survey system 



P48 survey characteristics 



Telescope 


Palomar 48-inch (1.2m) Samuel Oschin 


Camera field dimensions 


3.50x2.31 degrees 


Camera field of view 


8.07 square degrees 


Light sensitive area 


7.26 square degrees 


Plate scale 


1.01 arcsec / pixel 


Efficiency 


66% open- shutter (slew during readout) 


Sensitivity (median) 


m R ~20.6 in 60 s, 5<r 




m g /~21.3 in 60 s, 5<r 


Image quality 


2.0 arcsec FWHM in median seeing 


Filters 


g' & Mould-R; other bands available 


P48 survey camera CCD array 


Component CCDs 


12 CCDs; 1 non-functional 


CCD specs 


2Kx4K MIT/LL 3-edge butted CCDs 


Array Leveling 


Flat to within 20 microns 


Pixels 


15 microns/pixel; 100% filling factor 


Chip gaps 


Median 35 pixels (35 arcsec) 


Readout noise 


<20 e" 


Readout speed 


30 seconds, entire 100 MPix array 


Linearity 


better than 0.5% up to 60K ADUs 


Optical distortion 


maximum 1" at array corners 




compared to flat grid 


P60 followup camera specifications 


Telescope 


Palomar Observatory 60-inch (1.5m) 


CCD specs 


2Kx2K CCD 


Plate scale 


0.38 arcsec / pixel 


Readout noise 


5 e~ (amp 1); 8 e~ (amp 2) 


Readout speed 


25 seconds (full frame) 


Linearity 


better than 1% up to 20K ADUs 


Sensitivity (median) 


mj,~2 1 .6 in 120 s, 5<x 




nv ~21 .3 in 120 s, 5<x 




mi/«21.1 in 120 s, 5<r 




m z , -20.0 in 120 s, 5cr 



We summarize here the upgrades and provide on-sky per- 
formance test results for the PTF survey camera on P48; more 
detail on the camera engineering is provided in lRahmer et al.l 
(2008). TableCD summarizes the survey and P60 followup sys- 
tem specifications. 

2.1. CCD Array 

The core of the PTF survey camera is a 12K x 8K mosaic 
made up of twelve 2Kx4K MIT/LL CCID20 CCDs, arranged 
in a 6x2 array (Figure[3]). Three of the CCDs are high resistiv- 
ity bulk silicon (HiRho); the rest are standard epitaxial silicon 
(EPI). The EPI chips reach QEs of approximately 70% at 650 
nm; the HiRho devices have somewhat higher QE, reaching 
^90% at 650 nm (FigureSJ). The HiRho chips also have lower 
fringing levels due to their increased thickness. 

Depth of focus, and thus the flatness of the CCD array, is an 
important issue for the P48's f/2.5 beam. Approximately 92% 
of the PTF survey camera's CCD array is within 20 jim of the 
reference focal plane position. Taking a detailed error budget 
into account, including terms for telescope jitter, atmospheric 
turbulence and the optical quality, we predicted that 89% of 
the array would provide images meeting our 2.0"FWHM im- 
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Fig. 2. — Cut-away diagram of the PTF survey camera after modifications for use on the Palomar 48" telescope. The diagram is orientated so downwards is 
in the direction of the P48 primary mirror when mounted in the telescope. From bottom up, the components are: the filter exchanger, shutter, and the dewar 
assembly including the CCD window and CCD array. 
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Fig. 3. — Schematic of the PTF survey camera focal plane (with greatly 
expanded chip gaps). The HiRho chips (shaded) are high-resistivity silicon 
with improved red quantum efficiency (Figure [4). The read noise and gain is 
noted for each chip; some chips show relatively high read noise and are being 
optimized in ongoing work. CCD03 is unresponsive and will be replaced 
later in the PTF project. Chip gaps are noted in arcseconds, based on on-sky 
astrometric solutions. 
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Fig. 4. — The quantum efficiencies of the nine EPI chips in the CCD ar- 
ray (solid line) and the three HiRho chips (dashed line). Filter transmission 
curves for the two primary PTF survey filters are also shown (filled curves; 
green is SDSS-g' and re d is Mould-R). The QE curves and Mould-R trans- 
mission are adapted from Cuillandre et al. ( 2000). 



age quality budget and that the remainder of the array would 
be within 0.2" of the specification. On-sky tests confirm this 
( 32. 9i the best images so far taken have a median FWHM of 
1.5"across the array). The instrument tip and tilt was adjusted 
to better than 0.003 degrees tilt across the entire array using 
20 yum shims at the dewar mounting points and on-sky image 
tilt measurements from FWHM maps. The final full-field-of- 
view image quality is better than 2.0 arcsec FWHM in median 
seeing ( 32. 91) . 

The CCDs are 3 -edge buttable and are separated by gaps of 
^500 microns, corresponding to 35 arcseconds (see Figure [3] 
for more details). The CCD array itself has a low number of 
bad pixels and offers excellent image quality across the field 
of view (Figure [5} . 

The full CCD array readout time is 30.7 seconds, using 
independent amplifiers for each chip and two Generation 3 
Leach (SDSU) devices reading out 6 CCDs each. This read- 
out time is approximately half that of the original camera elec- 
tronics and was specified to match the expected telescope slew 
and settle time for PTF operations. The improved speed was 
reached at a small cost in readout noise; even so the final read- 



out noise is <15e" on all chips, well below the expected sky 
noise. 

During the camera upgrade process CCD03 became unre- 
sponsive. After diagnosing the fault as a problem inside the 
CCD package itself we decided that attempting to repair the 
CCD would be both an unacceptable risk to the array and an 
unacceptable delay to the project. We are continuing to evalu- 
ate ways to repair or replace the malfunctioning CCD, but for 
at least the first year of science operations PTF is accepting 
the 8% loss of imaging area. 

2.2. Optics 

We replaced the original flat dewar window with a field flat- 
tening optic designed for P48's curved focal plane. The new 
optics are designed to provide better than 2" image quality 
over the entire 3.4-degree span of the image in median see- 
ing, including atmospheric, optical and mechanical tracking 
error terms. The optical distortion is 0.11% at the corners of 
the array relative to a uniform grid, a level comparable to dif- 
ferential atmospheric refraction over our field of view (Figure 
[6] shows the distortion field). 
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Fig. 5. — A typical R-band image taken with the PTF survey camera on the Palomar 48" telescope, with the same orientation as Figure [3] The image has been 
bias-corrected and flat fielded but, to show the cosmetic quality of each chip, no dithering has been applied. 



The telescope optics were also optimized: the primary mir- 
ror supports were tuned; all optical surfaces were cleaned or 
refurbished; baffles were improved to control stray light; the 
telescope tube was made light tight; and internal stabilized 
calibration sources were added to allow daytime quality con- 
trol. 

2.3. Cooling system 

The original liquid nitrogen cooling system was replaced 
with a closed-cycle cooling system, implemented via a com- 
pact cold head located close to the CCDs. A Poly cold Com- 
pact Cooler (formerly known as CryoTiger) with an expansion 
head with no moving parts was chosen to minimize vibration. 
We attached an H- shaped copper head spreader to the cold 
head to provide a nearly isothermal surface to which to attach 
flexible copper heat exchange straps leading to the CCD array 
assembly. A heating system attached to the dewar window 
actively maintains the CCD array temperature at 175 K. 

2.4. Shutter 

The shutter is a commercial unit supplied by Scientific In- 
strument Technology, and employs dual split blades to avoid 
beam obscuration by the blades when the shutter is open. The 
blade transit time is 1 .300 seconds, but the blades can be inde- 
pendently controlled to allow shorter exposure times in a slit 
scanning across the detector. 

2.5. Filter changer 

The filter changer assembly is a custom-built mechanism 
that provides for motorized selection of one of two filters; one 
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Fig. 6. — Distortion map across the PTF survey camera focal plane, relative 
to a flat grid on the sky and based on Zemax modeling of the optical system. 
A 20" distortion scale is indicated at the bottom of the figure. The maximum 
distortion at the corners of the CCD array is 7", or 0.1 1%. 



filter is always in the field. Filter exchanges are completed 
in 15 seconds, using a micro-stepper motor. Absolute filter 
positioning is registered by running the filter frame against a 
mechanical hard stop. Normal PTF operation is expected to 
use at most two filters in each night, but the system has been 
designed for easy manual swapping of filters as necessary. 

2.6. Camera Electronics and Software 
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The PTF survey camera array of twelve CCDs is divided 
into two banks of six CCDs, with each bank handled sepa- 
rately by its own detector controller. The controllers are Gen- 
eration 3 Leach (SDSU) devices comprised of three 2-channel 
video boards, three clock boards, and a timing board. Com- 
munication between a controller and its respective host com- 
puter is conducted via a fiber optic serial link between the 
timing board and a corresponding PCI card in the host com- 
puter. The controllers share a common clock so that the ex- 
posures and readouts of the two controllers are synchronized 
with each other. The master camera computer also handles all 
the non-detector camera hardware - the shutter, filter changer, 
temperature sensors and calibration LEDs. The camera con- 
trol software is custom-built for the PTF survey camera, and 
is based on pan VIE W, the Pixel Acq uisition Node (PAN) ver- 
sion of Arc VIEW (lAshe et al.ll2002h . 

2.7. PTF Observatory Control System and Scheduler 

The survey operations are performed robotically by the P48 
Observatory Control System (OCS), written in MATLAB. 
The OCS is responsible for controlling the camera, filter ex- 
changer, shutter, telescope, dome, and focuser, on the basis of 
feedback information from the scheduling system, the camera 
a nd te lescope, a weather station, and the data quality monitor 

The OCS is responsible for sequencing all PTF observa- 
tions, starting with the bias frames before sunset, through fo- 
cus images, and finally the science images. The PTF sched- 
uler is responsible for selecting the next target for observa- 
tion. The PTF schedule is not pre-defined; the next target is 
selected 30-90 s before each exposure starts, during the previ- 
ous exposure. 

The criteria for selection of the next target include: (i) Sun 
altitude; (ii) excess in sky brightness due to the Moon, us- 
ing the algorithm of lKrisciunas & Schaeferl(ll99l[) : (iii) Moon 
phase; (iv) time needed to move telescope to the next target; 
(v) time needed to move dome to the next target; (vi) air- 
mass; and (vii) time since the last sequence of observations 
of this field. I n add ition for some programs, such as the 5- 
day-cadences ( 33.21) . the scheduler is required to optimize the 
field selection so that each field is observed twice during the 
night for asteroid detection. 

These criteria are used to calculate weights for each field, 
and the scheduler selects the target with the highest weight for 
observation. In order to keep the operations of the scheduler 
as simple as possible most of the weights follow a step func- 
tion. An exception is the cadence (time since last observation) 
weight (W) which is of the form: 



l+exp({[/Z)-/A]-T}/s) 

where JD is the current Julian day, JD\ is the Julian day cor- 
responding to the last successful sequence (pair of images) of 
observation, r is the cadence, and s is a softening parameter 
that defines the time scale in which the weight goes near to 
near 1 . 

2.8. Data Quality Monitoring 

A software suite running at the mountain, the Data Qual- 
ity Monitor (DQM), analyses all the recorded images as 
they are taken. An initial check searches for problems with 
file corruption, truncation, or missing FITS headers. Next, 
the image is searched for stellar objects [using SExtractor 



(Bertin & Arnouts 1996)], and the average FHWM and other 
image quality parameters are checked for obvious problems. 
Finally, if it has passed all the checks, the file is sent to the 
PTF pipelines. 

To detect more subtle image quality problems such as in- 
creased electronic noise the P48 twice-daily automatically 
takes test exposures of the telescope tube cover. A band of 
temperature- stabilized LEDs provides illumination stable to 
<0.1%. Because the illumination has passed through every 
optical element of the system, comparing the test exposures 
to reference exposures provides a simple test if anything in 
the telescope + camera + electronics systems is not operat- 
ing as expected. Many problems can be found and corrected 
before the start of nighttime operations. 

Finally, both the LNBL and IPAC pipelines produce image 
quality and sensitivity measures for all detected fields, pro- 
viding a long-term monitoring facility to allow us to optimize 
our observing setup and equipment. 

2.9. On- sky Performance 

The PTF survey camera system achieved first light on 13 
Dec 2008. Initial commissioning work was completed in Jan- 
uary 2009 and the system has been used for initial science 
test operations since then (with small gaps for engineering). 
Here we describe on- sky performance results from the first 
two months of operation. Since the system continued to be 
commissioned during this time, the performance is likely to 
improve as the PTF survey continues. 

A total of 91.7% of the pixels in the full 12Kx8K array are 
light-sensitive and not dark, occulted or otherwise bad (this 
increases to >99% good pixels if we exclude the unrespon- 
sive CCD). Most of the approximately 200 bad columns in 
the array are concentrated on CCD05 (at one corner of the 
field). CCD09 offers the best combination of PSF quality and 
CCD cosmetics. The chips are linear to better than 0.5% up 
to 60,000 ADUs (typical gain is 1.6e" ADIT 1 ). 

During the design of the PTF survey camera we specified 
a FWHM of 2.0" during median Palomar conditions (1.1" 
seeing). This requirement was based on a number of fac- 
tors - the minimum PSF size that is reasonably sampled by 
our 1 "/pixel plate scale; what could be reasonably achieved 
with the P48 optical setup; optimization of the survey limiting 
depth; and the requirement to detect supernovae in the cores 
of bright nearby galaxies. Analysis of the first two months 
of PTF on- sky data shows that we have met or exceeded this 
requirement. Typical images taken in good seeing show both 
sub-2.0" FWHMs and very small variations in image quality 
across the entire field of view (Figures [7] & [8]). 

The survey depth is correspondingly encouraging (as shown 
in Figure O. Based on all the data taken during the last 
two months of commissioning, the system provides median 
5<x, 60-sec limiting dark-time magnitudes of m g >=21.3 (25 th 
percentile m g /=21.7; 75 th percentile m g /=21.0) and niR=20.6 
(25 th percentile m R =20.8; 75 th percentile m R =20.3). 

The entire PTF survey system (the P48 telescope, PTF sur- 
vey camera, OCS and the DQM) has a current uptime of 
^90%, excluding losses due to weather. As initial science 
operations proceed, many issues are being addressed and the 
uptime is continually improving. The survey efficiency is cur- 
rently approximately 50% open-shutter time, compared to a 
goal (and theoretical maximum) of 66%. This is due to the 
cumulative effect of many small inefficiencies in the system 
operation and is continually improving. 
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Fig. 7. — A map of the FWHM variations across the PTF survey camera fo- 
cal plane array from a typical R-band observation (with the same orientation 
as Figure 0. Image FWHM measurements are interpolated over the unre- 
sponsive CCD for ease of interpretation. The residual FWHM variations are 
on the order of 0.2" and are due to the 20//m flatness of the CCD array and 
slightly decreasing optical performance at the edge of the field of view. 




Fig. 8. — 100" cutouts from R-band PTF survey camera images, from the 
first set of PTF optical transient detections. Crosshairs show the transient 
detection position, most inside galaxies. The properties of the detected tran- 
sients are given in §6\ with the exception of PTF-OT3, which is a re-discovery 
of SN2009an. The images are distributed widely across the focal plane. 



3. THE PTF OBSERVING STRATEGY 

PTF will pursue four distinct experiments optimized for 
different times of the year. Each project is designed to use 
P48 only as the survey telescope while other PTF resources 
are reserved for followup observations. We here describe the 
observing strategy for each of the major PTF experiments; 
iRau et al.1 (|2009[) give a detailed discussion of the science 
planned for each project. The experiments are summarized 
in Table [2 
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Fig. 9. — Completeness limits from a typical g'-band, photometric, PTF 
survey night. Each curve represents the minimum detectable magnitude at 
either 50% or 95% completeness; each exposure of each chip is treated as an 
individual image. The spread in values arises from variations due to seeing, 
the imaging quality of each chip and the airmass of the observations, and thus 
represents typical performance for the actual PTF observation strategy. 



TABLE 2 

The PTF experiments planned for the first year of operations. 



Experiment 


% of total 


Cadences 


Filter 


Months 


5 -day cadence (5DC) 


41 


5d 


R 


Mar-Oct 


Dynamic cadence (DyC) 


40 


1 min - 5 d 




all year 


Orion 


11 


1 min 


R 


Nov-Jan 


Ha 


8 




Ha 


all year 



3.1. Observing Setup & Filter Selection 

The three major PTF experiments, together accounting for 
92% of the project's P48 time, are based on two filters: SDSS 
g' and Mould-R (Figure |4]). These filters were chosen to op- 
timize the survey detection limits in full moon (R-band) and 
dark sky (g'-band) conditions for a variety of different tran- 
sient classes (see IRau et al.ll2009h . During dark-sky condi- 
tions, the reduced quantum efficiency in the g' filter is more 
than compensated for by the darker sky and the blue colors of 
many of our targeted transients. Mould-R is very similar to 
SDSS r' and allowed us to reuse an existing CFH12K filter; it 
is the standard filter for programs that are expected to spend 
at least part of their time operating while significant moon- 
light is present. Because no followup is done with P48, we 
elected to use standard passbands for the survey to allow a 
meaningful comparison of PTF detection magnitudes to fol- 
lowup telescope results. 

The initial standard exposure time for the 5DC and DyC 
surveys is 60 seconds, based on the readout time of the cam- 
era, expected system efficiency, and the expected sky noise 
and dark current in the exposures. In later stages of the PTF 
project this exposure time will be optimized for the different 
programs expected to be pursued. 

3.2. 5 Day Cadence 

The 5 -day cadence experiment (5DC) will run yearly from 
March 2nd until October 30th using 65% of the time in that 
period. The main goal of this experiment is to construct large 
samples of type-la SNe and core-collapse SNe. It will also 
study AGNs, Quasars, Blazars, CVs, extragalactic novae, lu- 
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minous red novae, and tip of the AGB variability. The exper- 
iment will observe a footprint of approximately 8000 square 
degrees, with \b\ > 30°, ecliptic latitude \j3\ > 10°, and with a 
median cadence of 5 days. At a given time, the active area for 
the 5DC search (including all overheads) will be 2700 square 
degrees. The 5DC observations will be conducted in almost 
all lunar phases in R-band. In each epoch we will obtain two 
60 s exposures separated by 60 min. The two images will 
be used to remove cosmetic artifacts, cosmic rays and to find 
solar system objects. 

3.3. Dynamic Cadence 

The dynamic cadence experiment (DyC) is designed to ex- 
plore transient phenomena on time scales shorter than ~5 days 
and longer than one minute. The experiment will explore un- 
known territories in the transients phase space - mainly fast 
transients with time scales of less than a day. Therefore, the 
dynamic cadence will be an evolving experiment. Every six 
months the PTF collaboration will review the results of the 
dynamic cadence experiment and decide how to continue. 

Motivated by a search for short-lived transients in the 
brightness gap between novae and supernovae(-10<M< 
-16), the DyC experiment is targeted at fields that are lumi- 
nosity concentrations in the local universe. Since these tran- 
sients are expected to be fainter than M 16, the search vol- 
ume is limited to DM (distance modulus) < 36.5 - i.e. less 
than 200 Mpc. 

The total luminosity of nearby galaxies in the PTF FoV was 
computed on a fine grid of possible pointings (limited by dec- 
lination > -30). This was done seven times with seven differ- 
ent maximum distances (in steps of 0.5 between 30.5 < DM 
< 36.5). We assigned a weight for each pointing based on the 
fraction of maximum luminosity in that list. A total weight 
for each pointing was computed by averaging the individual 
weights. Finally, pointings with the highest weights were cho- 
sen for PTF. This list naturally includes the nearest, brightest 
galaxies (e.g. M31, M51, M81) and nearest galaxy clusters 
(e.g. Virgo, Perseus, Coma), but also covers the other acces- 
sible nearby mass concentrations from physical and chance 
galaxy alignments. 

This experiment will detect RR Lyr stars in the Galactic 
halo, CVs, AM CVn stars, flare stars, SNe, eclipsing binaries, 
and solar system objects. Moreover, this experiment will put 
strong constraints on the existence of any optical counterparts 
that may be associated with radio transients. 

3.4. Orion Field 

The Orion project has been assigned 40 consecutive nights 
per year for three years to perform intensive time- series ob- 
servations of a single field in the Orion star forming region, 
with the aim of detecting close-in Jupiter-sized planets transit- 
ing young stars. The observations will be acquired in R-band 
with no dithering in a continuous sequence of 30- second im- 
ages, with the goal of producing high precision (better than 
1 %) differential photometry. The field will be chosen to rep- 
resent a stellar age range where protostellar discs should be on 
the point of dissipation (5-10 Myr) - giving important clues 
about planet formation and migration - and to optimize for 
source density and avoid regions where the stellar variability 
in the ensemble is severe enough to prohibit precise differ- 
ential photometry. The IPAC PTF pipeline will be used for 
standard image reduction, source identification, and astrom- 
etry, but a dedicated photometric pipeline is being developed 
to produce the high-precision differential photometry. 



3.5. Deep Ha Sky Survey 

For three nights each lunation during full Moon we will 
carry out a 3n sr sky survey in Ha. The survey area will in- 
clude all the sky with S > -25°. Each field will be observed 
twice in four narrow-band (approximately 10 nm width) filters 
around the Ha line covering the redshift range of about zero 
to 0.05, corresponding to distances of to 200 Mpc. Narrow- 
band filters that give a consistent bandpass across the entire 
field of view of the camera are difficult to produce, especially 
in the fast beam of the P48. However, a description of us- 
able Ha filters fo r a Schmidt telescope with the sam e optical 
design is given in Parker & Bland-Hawthorn (1998), and our 
industrial partners are confident that these filters can be con- 
structed. 

We estimate that the Ha survey sensitivity limit will be 
^0.6R in a 2.5 arcsec aperture 1 ; a detailed com parison of 
this se nsitivity to other Ha surveys can be found in llRau et all 
(2009). We are currently designing a scheme that will allow 
us to photometrically calibrate these data to about 10% accu- 
racy. 

4. DATA REDUCTION AND TRANSIENT DETECTION 

All PTF data taken at P48 is automatically routed to 
two pipelines: a realtime transient detection pipeline, opti- 
mized for rapid detection of interesting objects, and a longer- 
term archival pipeline optimized for high precision and easy 
searching of the results. 

4.1. Realtime Transient Detection 

The PTF quick-reduction and subtraction pipeline was de- 
signed from the start to take full advantage of the paral- 
lel high performance computing facilities at the National 
Energy Research Scientific Computing Center (NERSC, 
http://www.nersc.gov). The major improvement to the hard- 
ware used in this pipeline over previo us versions designed by 
the Supernovae Cosmology Project (Perlmutter et al. 1999) 
and Nearby Supernovae Factory (lAldering et al.l |2002|) at 
LBNL is the use of the NERSC Global Filesystem (NGF). 
NGF is a 300TB shared filesystem with over 1 GB/sec band- 
width for streaming I/O which can be seen by all of the high 
performance computers. On the software front, the improve- 
ments include a tight coupling of a PostgreSQL database in- 
volved in tracking every facet of the image processing, ref- 
erence building, image subtraction and candidate detection 
along with a complete re- write of the processing and subtrac- 
tion codes, described below. 

Incoming P48 images are immediately backed up on the 
High Performance Storage System, a 3 PB tape archive sys- 
tem at NERSC. Processing begins by splitting the packed 
images apart by chip, applying crosstalk corrections, and 
performing standard bias/overscan subtraction and flatfield- 
ing. After this point, all operations are performed on a 
chip-by-chip basis in parallel. Catalogs of sources are 
created for each image via the Terapix SExtractor code 
(Bertin & Arnouts 1996), which then is fed to the astrome- 
try.net (http://astrometry.net) code to perform an astrometric 
solution. At this point a comparison to the USNO catalog is 
made to determine the zeropoint and 3-cr limiting magnitude 
of the image along with a calculation of the seeing. The im- 

1 R is a unit of surface brightness commonly used in aeronomy. One R (one 
Rayleigh) is 106/(47r) photons per square centimeter per steradian per second. 
For the Ha line the intensity in cgs units is 2.41 x 10~ 7 erg cm -2 s _1 sr -1 . 
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Fig. 10. — The detection images of SN 2009av, PTF's first confirmed tran- 
sient detection. The top and bottom row show the same images at different 
stretches. 



age is then loaded into the processed image database, storing 
all relevant information from the FITS headers. 

After several images are obtained for a given PTF pointing, 
we are able to create a reference image. The reference im- 
ages are crea ted via the Terapix software Scamp and Swarp 
(Bertin 2006) after querying the processed image database in 
order to obtain the highest quality input images for a given 
pointing/chip combination. The references are assigned ver- 
sion numbers and (along with their relevant information) are 
stored in the reference database. 

Subsequent new images are processed as above, and if a 
reference image exists for them, a subtraction is performed. 
First, the new image is astrometrically aligned to the refer- 
ence by utilizing Scamp and the new and reference catalogs. 
The reference is then Swarped to the size and scale of the new 
image and a subtraction is performed using HOTPAnTS 2 . The 
subtraction is performed in both the positive and negative di- 
rection (reference minus new and new minus reference) to de- 
tect both positive and negative flux changes. Candidate tran- 
sient sources are detected via SExtractor and along with other 
relevant parameters (location with-respect-to bad pixels, etc.) 
are stored in a candidate and subtraction database. The sub- 
traction image that lead to the first PTF supernova discovery 
is shown in Figure [TOl 

In addition, PTF takes advantage of DeepSky 
(http://www.deepskyproject.org). DeepSky was started 
in response to the needs of several astrophysics projects 
hosted at NERSC. It is a repository of digital images taken 
with the point-and-stare observations by the Palomar- 
QUEST Consortium a nd the Near Earth Asteroid Team 
(Djor govski et al.l 120081) . This data spans nine years and 
more than 15,000 square degrees, with 20-200 pointings on 
a particular part of the sky and cadences from minutes to 
years. For a large fraction of the survey DeepSky achieves 
depths of mR > 23 magnitude. In total there are more 
than 11 million images in DeepSky. This historical data 
set compliments the Palomar Transient Factory survey by 
identifying known variable stars, AGN, and the detection of 
low surface brightness host galaxies of supernovae. 

During commissioning human scanners were utilized to re- 
ject "junk" transient detections via a web-interface. The sub- 

2 http : //www . astro . Washington . edu/users/becker/hotpants . html 



traction system and rejection algorithms are currently being 
improved with an aim of full automation by making use of 
machine learning techniques. The results of human scanners 
are being used to train an automatic classifier; a more detailed 
account of this classification technique is in preparation. 

Any statistical inference of PTF transient rates will require 
a careful measurement of the survey detection efficiency. We 
have constructed a pipeline that produces realistic fake vari- 
able objects and adds them into PTF images, on which we 
apply the standard PTF search algorithms to determine the 
detection efficiency for any specific data subset. We cur- 
rently model three generic classes of transients: variable stars, 
SN/AGN-like sources (i.e., point sources in host galaxies), 
and "rogue" transients (point sources on random positions in 
the image). The transients are directly injected into PTF im- 
ages using PSFs modeled from nearby stars in the fields. 

4.2. Realtime Automated Transient Classification 

The Transients Classifi cation Pipeline (TCP; iBloom et~ai1 
120081: IStarr et al.l l2008allbh is a parallelized, Python-based 
framework created to identify and classify transient sources 
in the realtime PTF differencing pipeline (Section [47Tb . The 
TCP polls the candidate database from that pipeline and re- 
trieves all available metadata about recently extracted sources. 
Using the locations and uncertainties of the transient candi- 
date objects, the TCP either associates an object with ex- 
isting known sources in the TCP database or, after pass- 
ing several filters to exclude non- stellar (e.g. known mi- 
nor planets) and non-astrophysical events, generates a new 
so urce. The Bayesi an framework for event clustering is given 
in Bl oom et al] (120081) . An other classification e ngine will also 
be used on PTF transients (Maha bal et al .112008 alibi) . 

Once a transient source has been identified, the TCP gener- 
ates "features" which map contextual and time-domain prop- 
erties to a large dimensional real-number space. After gen- 
erating a set of features for a transient source candidate, the 
TCP then applies several science classification tools to deter- 
mine the most likely science class of that source. For rapid- 
response transient science, a subset of features — such as 
those related to rise times and the distance to nearby galax- 
ies — are most useful. As light curves are better sampled and 
colors are obtained more features are used in the classifica- 
tion. The resulting science class probabilities are stored in a 
database for further data mining applications. Sources with 
high probabilities of belonging to a science class of interest to 
the PTF group will be broadcast to the PTF's "Folio wup Mar- 
shal" (Section[5]) for scheduling of followup observations. As 
more observations are made for a known source, the TCP will 
autonomously regenerate features and potentially reclassify 
that source. 

4.3. The IPAC PTF Image Processing Pipeline and Data 
Archive 

While the transient pipeline will be optimized for rapid 
turn-around, the IPAC pipeline will be geared towards pro- 
viding the best possible photometry consistent with hardware 
capabilities and data rate. To enable the generation of optimal 
calibration files (primarily flat fields), the IPAC pipeline will 
commence image processing only after all of a given night's 
observations have been received. Additional source associ- 
ation operations, database updates, and data transfer to the 
Infrared Science Archive (IRS A (http://irsa.ipac.caltech.edu) 
will be carried out over a period of several days. Processed 
images and incremental source lists will be available through 
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IRSA within 2 days from the time the observations were 
taken. 

After a given night's observations have been completed 
and received at IPAC (as indicated by receipt of end-of-night 
image manifest), the image processing pipeline is initiated. 
Parallel processing on 11 multi-processor Linux machines 
is employed to simultaneously process the PTF mosaic's 11 
working detectors and meet the data rate requirement. Key 
pipeline-processes are multi-threaded to take full advantage 
of the 8 -core CPU in each of the Linux machines. 

The pipeline initially generates night- specific bias frames, 
and trims and bias- subtracts all images in the usual manner. 
Nightly flat fields are generated for each detector by source- 
masking and median-combining all of the night's sky observa- 
tions. This is an involved and memory intensive process and, 
given the number of exposures involved, takes 3-4 hours with 
the processing hardware currently in place. On completion, 
these flat fields are used to process the entire night's science 
observations. Raw and processed images as well as all ancil- 
lary and calibration frames are subsequently passed to IRSA 
for archiving and retrieval on demand. 

Source detection is carried out using SExtractor 
(Bert in & Arnoutsl Il996|) . SExtractor is applied to each 
image a number of times in the pipeline to mask sources prior 
to flat field generation, measure the effective seeing, generate 
pixel weights and masks, and finally to generate photometry 
for all detected sources. Pixel masks are generated to flag bad 
rows and columns, cosmic rays, saturated pixels, and charge 
bleeding. The final application of SExtractor is carried out 
after multiplying the processed image by a pixel area map to 
preserve point source photometry. 

For each image, the SExtractor output list is then astr omet- 
rically calibrated using the Terapix Scamp software (iBertinl 
2006). Photometric zeropoints are determined b y calibrat- 
ing ag ainst stars in the Sloan Digital Sky Survey (lYork et al.l 
2000), taking color terms into account. With appropriate air- 
mass corrections, these zeropoints are then be applied to all 
sources outside the SDSS footprint. Because the PTF sur- 
vey camera's CCD conserve signal when stars are bloomed 
due to saturation, absolute photometric calibration may also 
be carried out using bright (my~l 1) sa turated Tych o stars for 
which the griz magnitudes are known ( Ofek 2008). The den- 
sity of such stars in the PTF survey camera field of view is 
high enough and we estimate it will be possible to achieve 
accuracies of better than 5% in both the g' and R bands in 
non-SDSS fields. 

Currently in development is a source association pipeline 
that will match, on a nightly basis, each detected source with 
all previous detections of that source. Appropriate pointers 
will be incorporated into the IRSA database, enabling users 
to query the photometric time history of every object ever de- 
tected by the PTF." 

5. FOLLOWING-UP DETECTED TRANSIENTS 

The automated Palomar 60-inch is the primary workhorse 
for photometric follow-up of PTF transients. Further photo- 
metric and spectroscopic characterization is performed on a 
worldwide network of other telescopes, including the Palo- 
mar Hale Telescope, the Keck telescopes, PARITEL, and the 
LCOGT telescope network. 

5.1. Photometric Followup with P60 

The automated P60 observation queue dCenko et al.| [2006) 
is sent all candidates flagged as bona fide astrophysical tran- 



sients (i.e. a new point source detected at least twice in P48 
imaging and which either lacks a counterpart or has a galaxy 
host in the reference imaging). A 120 s snapshot in g',r',i', z' 
is obtained for quick classification. On a typical night with 
median weather conditions, this reaches a limiting magni- 
tude of m^ *21.4, m; *21.2, mf ^20.9 and m' z «19.8. Our 
efficiency currently allows twenty two-minute exposures per 
hour. 

As soon as the 60-inch data is obtained, an automated 
pipeline detrends the data and solves for an astrometric and 
zeropoint solution. Next, if the location is in SDSS, the 
pipeline downloads the corresponding SDSS image (mo- 
saiced to fit the P60 FoV) to serve as a template. A convo- 
lution kernel is measured and image subtraction performed. 
On the subtracted image, PSF photometry is used to measure 
the magnitude at the source location and the result is sent to 
the PTF database. If the source is not detected, an upper limit 
is returned. Currently, if the location is outside SDSS, direct 
photometry is performed. We note that host galaxy light can 
be a significant contaminant which is not properly corrected 
for in case of direct photometry. Hence, even for fields out- 
side SDSS, we are investigating other options for template 
imaging for image subtraction. 

The 60-inch snapshot observation provides independent 
confirmation that the candidate is a real transient (and not a 
subtraction artifact), a color for quick classification (for ex- 
ample, m g - m r « suggests that it is a SN la), and a second- 
epoch magnitude to assess photometric evolution. All this 
information helps direct further spectroscopic and multi-band 
photometric follow-up with P60 and other telescopes. 

5.2. The Followup Marshal 

We have designed a centralized follow-up Marshal to en- 
sure that our discoveries are efficiently distributed to our 
follow-up resources and that the temporal coverage of a given 
target meets the requirements of the appropriate science pro- 
gram. The design of this system is complete, but the auto- 
mated followup program has not yet begun. A later paper will 
detail the performance of this novel system. 

The first duty of the follow-up Marshal will be to moni- 
tor the output of the transient classifier to identify any new 
sources detected by P48 that warrant follow-up observations. 
For this purpose, the Marshal will maintain a list of target de- 
scriptions (in terms of the classifier output) for each science 
program to define trigger conditions for same night and long 
term monitoring. For example, a supernova program may de- 
fine same night triggers for all objects with at least a 25% 
chance of being a supernova but no more than a 1 % chance of 
being an asteroid. Discoveries worthy of same night follow- 
up will be sent primarily to P60, although sources identified 
at the end of the night may fall to Faulkes Telescope North 
(FTN; Haleakala, Hawaii) in rare cases where twilight pre- 
empts the required observation. These requests will be issued 
as ToO triggers, similar to the gamma ray burst triggers al- 
ready handled by both P60 and FTN except that interruption 
of integrations or read outs will not be required. 

After completing an observation, the data must be analyzed 
to extract the target parameters (e.g. the magnitude and posi- 
tion in photometric data or the spectral classification, redshift, 
etc. for spectral observations). The responsibility for this data 
reduction lies with the partner in charge of the instrument, and 
the high level data will be sent back to the Marshal. The Mar- 
shal will then relay this information back to the classifier to 
refine the classification, and it will archive the observation in 
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Fig. 11. — Keck LRIS spectra of several different types of optical transients discovered and classified by PTE Spectra are not redshift-corrected, and are 
vertically offset for clarity. Details of the detected transients can be found in Table[3] 



its own database for future use. 

6. THE FIRST PTF TRANSIENT DETECTIONS 

We conclude with the first results from the PTF sur- 
vey. PTF's first bona fide transient was detected 
at ny = 18.62±0.04 mag on UT 2009 March 2.347 
( Kul karni et al.l l2009|) after subtraction of a reference tem- 
plate made from PTF images acquired on Feb 25, 26, and 
28 (Figure [TO]). The transient, at RA = 14:23:55.82 Dec = 
+35:11:05.2 (J2000, uncertainty <1"), was offset 0.9" W and 
2.7" S from the core of an SDSS galaxy with m g = 17.78 and a 
redshiftof z=0.0555. 

We obtained spectra of the transient on March 8.388 UT 
with the Double Beam Spectrograph (DBSP) instrument on 
the 5. 1-m Palomar Hale telescope. The transient turned out to 
be a normal Type la supernova at a redshift consistent with the 
apparent SDSS host. In this restframe, the expansion velocity 
derived from the minimum of the Sill (rest 635.5 nm) line was 
about 10,700 kms -1 . The best match found with t he Superfit 
SN spectral identification code (iHowell et al .120051) was to SN 
1994D at 8 days prior to maximum light. 

During the commissioning process many further transients 
have been detected and followed up (Table [3j Figure fTTT) . 
Eleven were detected during an R-band search of 498 square 
degrees (iKasliwal et al.l l20Q9bb . using reference images ac- 
quired by PTF prior to March 12. All eleven transients were 
spectroscopically classified on UT 2009 Mar 20 using DBSP. 
Twenty-five optical transients were found in a search of sev- 
eral nights of commissioning data taken between March 17 
and March 28, photometrically followed up with P60, and 
spectroscopically c lassified using the H ale, Hobby-Eberly and 
Keck telescopes dOuimby et al.ll2009al). A very bright super- 
nova was detected on April 17 dOuimby et alj|2009bl) and a 
further 15 opti cal transients were di scovered and followed up 
in May 2009 (IKasliwal et al.ll2009ab . These detections repre- 
sent a small fraction of the ultimate PTF discovery power, as 
the majority of the commissioning time was used for acquir- 



ing reference images. 

These initial transient detections have opened the Palomar 
Transient Factory for business. The full PTF survey is planned 
to start in June 2009 and the project is scheduled to continue 
until at least the end of 2012. 
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TABLE 3 

Optical Transients Discovered during PTF commissioning 



Name RA(J2000) DEC(J2000) Date UT 2009 Mag z 



Offset from Host Spect. Class. 



PTF09a 

PTF09d 

PTF09e 

PTF09f 

PTF09g 

PTF09h 

PTF09k 

PTF09o 

PTF09q 

PTF09r 

PTF09s 

PTF09t 

PTF09u 

PTF09v 

PTF09x 

PTF09y 

PTF09z 

PTF09aa 

PTF09ab 

PTF09ac 

PTF09ad 

PTF09aj 

PTF09as 

PTF09av 

PTF09aw 

PTF09bc 

PTF09bd 

PTF09be 

PTF09bh 

PTF09bi 

PTF09bj 

PTF09bw 

PTF09bx 

PTF09by 

PTF09ct 

PTF09cu 

PTF09dh 

PTF09ds 

PTF09ec 

PTF09fb 

PTF09fr 

PTF09fs 

PTF09fu 

PTF09go 

PTF09do 

PTF09ge 

PTF09gk 

PTF09gm 

PTF09gn 

PTF09ib 

PTF09ij 

PTF09jw 



14:23:55.82 
09 15 36.57 
09 14 34.16 

11 41 54.05 
15 16 31.50 

08 00 47.28 

12 26 17.84 
04 05 02.73 
12 24 50.17 
14 18 58.68 

12 06 51.67 
14 15 43.29 
14 29 39.61 

13 27 10.50 

13 21 45.15 
11 16 34.27 
11 54 42.23 

11 33 20.71 

09 22 15.69 

12 24 35.31 

11 03 06.64 
09 45 30.46 

12 59 15.85 

09 15 12.73 

14 15 19.36 

10 51 08.55 
08 07 29.72 

14 10 18.54 

12 24 39.20 

11 46 50.12 
11 18 06.46 

15 05 01.97 
14 30 50.42 

13 29 12.64 

11 42 13.88 

13 15 23.15 

14 44 42.08 
14 09 16.65 

13 12 54.49 

16 46 42.82 

14 50 00.12 

17 36 44.28 
16 33 10.88 

16 47 34.80 

17 38 26.04 

14 57 03.10 

15 06 11.06 
15 27 48.59 
15 29 10.95 
15 44 38.77 
14 32 14.62 

12 54 23.74 



+35:11:05.2 
+50 34 51.5 
+51 10 29.7 
+ 10 25 46.1 
+54 27 35.4 
+46 56 53.8 
+48 26 49.5 
+73 24 54.2 
+08 25 57.6 
+35 23 16.1 
+26 57 36.7 
+ 16 11 59.1 
+39 09 32.8 
+31 30 32.5 
+42 33 06.2 
+03 32 02.8 
+55 18 10.7 
-09 24 40.3 
+45 44 53.4 
+47 14 16.8 
+50 09 36.3 
+06 32 25.0 
+27 16 41.3 
+ 19 05 46.3 
+ 16 25 14.0 
+74 05 23.2 
+ 15 34 41.8 
+ 16 53 38.8 
+08 55 59.2 
+ 11 47 55.3 
+ 12 53 43.1 
+48 40 03.9 
+35 37 31.4 
+46 43 27.5 
+ 10 38 54.0 
+46 25 09.4 
+49 43 44.9 
+53 06 13.2 
+43 28 36.0 
+75 15 28.6 
+44 55 05.8 
+53 40 12.3 
+53 05 30.4 
+49 50 00.4 
+53 23 00.4 
+49 36 40.8 
+53 17 42.9 
+41 35 34.0 
+40 47 39.0 
+45 47 51.0 
+54 51 19.7 
+56 43 57.5 



Mar 02.3470 
Mar 17.1476 
Mar 17.1476 
Mar 17.1759 
Mar 17.2681 
Mar 17.1418 
Mar 17.1880 
Mar 17.1343 
Mar 17.2195 
Mar 17.2572 
Mar 17.1789 
Mar 17.2664 
Mar 17.2557 
Mar 17.2091 
Mar 21.3851 
Mar 21.1860 
Mar 21.4275 
Mar 21.1961 
Mar 21.3599 
Mar 21.4350 
Mar 21.4313 
Mar 21.2155 
Mar 25.1642 
Mar 25.1522 
Mar 25.2625 
Mar 20.2834 
Mar 20.1407 
Mar 26.2304 
Mar 28.4155 
Mar 25.1613 
Mar 25.1627 
Mar 28.4039 
Mar 26.2086 
Mar 26.2509 
Mar 27.1606 
Mar 26.2509 
Apr 17.375 
May 15.236 
May 15.259 
May 16.198 
May 17.172 
May 17.244 
May 17.247 
May 17.237 
May 15.220 
May 17.229 
May 17.253 
May 17.196 
May 17.196 
May 20.401 
May 20.288 
May 20.245 



18.6 

19.5 

20.0 

20.1 

18.4 

20.2 

20.9 

18.5 

19.6 

19.2 

17.8 

18.6 

20.0 

19.2 

20.2 

20.1 

20.1 

19.0 

19.5 

19.2 

19.8 

17.8 

19.1 

20.2 

19.7 

20.51 

16.9 

19.0 

20.5 

19.6 

18.1 

20.3 

18.8 

19.3 

20.3 

17.9 

20.3 

20.6 

19.8 

17.7 

19.1 

20.3 

20.5 

18.1 

20.2 

19.4 

20.0 

19.4 

20.4 

20.4 

20.3 

20.4 



0.06 0.9" W 2.7" S 



0.15 
0.15 
0.04 
0.12 
0.19 

0.09 
0.03 
0.05 
0.04 
0.13 
0.12 
0.25 
0.19 
0.19 
0.12 
0.17 
0.16 
0.20 
0.09 
0.19 
0.22 
0.17 
0.18 



0.18 
0.11 
0.14 
0.15 

0.10 

0.15 

0.06 

0.1 

0.179 

0.091 

0.042 

0.08 

0.109 

0.17 

0.047 

0.080 

0.064 

0.0 

0.082 

0.139 

0.123 

0.123 

0.153 



0.1"E 
0.5"E 
0.5"E 
0.4"E 
5.3"W 

1.8"E 

O.T'E 

0.2"W 

4.6"E 

0.0"W 

O.T'W 

0.8"W 

0.2"W 

0.1"W 

0.0"W 

0.0"E 

6.1"W 

0.6"W 

0. 2"W 

1. T'W 
2.0"W 
4.7"W 
3.9"E 

4.7"E 

1.6"W 

2.6"W 

0.1"E 

0.8"W 

4.2"E 

0.2"E 

5.6"W 

7.5"W 

0.9"W 

0.3"E 

0.0"E 

2.5"E 

1.4"W 

8.1"W 

0.2"W 

4.0"E 

0.1"E 

0.2"E 

1.1"W 

O.F'E 

3.7"W 

0.2"E 



0.8"S 
2.6"S 
4.5"N 
0.4"S 
6.3"N 

3.9"S 
0.8"N 
0.0"N 
1.4"S 
2.0"N 
0.5"N 
1.2"N 
2.1"N 
0.4"N 
0.2"N 
0.1"N 
1.2"N 
1.4"N 
1.5"N 
5.5"N 
2.9"S 
4.5"S 
1.9"N 

3.9"N 
0.8"N 
O.T'S 
1.3"N 
O.T'N 
3.3"S 
0.9"N 
5.2"N 
3.9"S 

0.2"N 
0.0"N 
0.2"S 
5.6"N 
3.5"S 
2.5"N 
0.2"S 
1.1"N 
0.1"S 

1.6"S 
0.4"S 
0.2"S 
2.1"N 
0.3"N 



SN la 
CV 
SN la 
SN la 
SN II 
SN la 
SN la 
CV? 
SN Ic 
SN II 
SN la 
SN II 
SN la 
SN la 
SN la 
SN? 
SN la 
SN la 
SN la 
SN la 
SN la 
SN la 
SN la 
SN la 
SN la 
SN la 
CV 
SN II 
SN la 
SN II? 
SN la 
SN Iln? 
SN? 
SN la 
SN II? 
SN II 
SN 
SN la 
SN la 
SN la 
SN la 
SN 
SN la 
SN II 
SN la 
SN? 
CV 
SN la 
SN la 
SN la 
SN 
SN la 
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